The formation and evolution of supermassive 
black holes and their host galaxies * 



Martin G. Haehnelt and Guinevere KauiTmann 

Max-Planck-Institut fur Astrophysik, Karl-Schwarzschild-Str. 1, 85740 Garching, 
Germany 



Abstract. We discuss constraints on the assembly history of supermassive black 
holes from the observed remnant black holes in nearby galaxies and from the emis- 
sion caused by accretion onto these black holes. We also summarize the results of 
a specific model for the evolution of galaxies and their central black holes which 
traces their hierachical build-up in CDM-like cosmogonies. The model assumes (i) 
that black holes, ellipticals and starburts form during major mergers of galaxies 
(ii) that the gas fraction in galaxies decreases with decreasing redshift (iii) that the 
optical bright phase of a QSO lasts for about 10^ years. The model succesfuUy re- 
produces the evolution of cold gas as traced by damped Lya systems, the evolution 
of optically bright QSOs, the remnant black hole mass distribution and the host 
galaxy luminosities of QSOs. 



1 Introduction 

The evidence for the existence of supermassive black holes has been steadily 
increasing over the last years. For the two most convincing cases our own 
galactic centre and NGC4258 ||l|,HJ^, the evidence is now beyond resonable 
doubt. The evidence that most nearby galaxies contain supermassive black 
holes is also compelling. Early suggestions of a linear relation between the 
black hole mass and the bulge mass have been corroborated by larger samples 
[Hj^. It is generally believed that we observe a significant fraction if not all 
the material falling into supermassive black holes by the radiation emitted by 
active galactic nuclei. A supermassive black hole therefore seems to "know" 
in which galaxy it will end up at the present day. This and the fact that su- 
permassive black holes contain as much as 0.2 to 0.6 percent of the baryonic 
mass of the galaxy , suggests that the the formation of stars in the bulges 
of galaxies and the assembly of supermassive black holes at their centre are 
closely linked . On the other hand, there is strong evidence that structures 
in the Universe form hierarchically, i.e. larger structures build up by merging 
of smaller structures. This is a generic feature of a wide class of structure for- 
mation scenarios, the so called cold dark matter (CDM) cosmogonies. Both 
galaxies and AGN activivity have been succesfully modelled within such hier- 
archical cosmogonies P,0,pdltl 2|Jl 3|,Q4,15l . Here we first review observational 
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constraints on the accretion history of supermassive black holes and discuss 
some clues for the formation mechanism. We then summarize the results of a 
specific model that describes the joint evolution of galaxies and supermassive 
black holes within a hierarchical cosmogony (see Kauffmann & Haehnelt 1998 
for more details ||l^). We assume here J7mat = 0.3, = 0.7, h — 0.65 and 
(78 = 1- 



2 The formation and evolution of supermassive black 
holes 

2.1 Black hole mass densities 

We can get some information on the assembly history of supermassive black 
holes in nearby galaxies by comparing the mass density of remnant black 
holes to that required to produce the radiation emitted during the accretion 
process [0. The mass density in remnant black holes can be inferred from 
the total mass density in bulges and the mass ratio of black hole mass to 
bulge mass. 
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The total mass density in bulges has been estimated by Fugukita Hogan and 
Peebles to be O.OOl/i^^ < l^buigo < O.OOS/i^^ |l|l. The normahzation of the 
bulge to black hole mass correlation is still a matter of debate. Magorrian et 
al claim a value of 0.6 percent while van der Marel argues that a value of 0.2 
percent is more realistic [^|J6|. We will adopt the latter value in the rest of 
the paper. 

The integrated emission by optically bright QSOs due to accretion onto 
supermassive black holes can also be used to infer the correponding mass 
density in supermassive black holes if an efficiency for the transformation of 
accreted rest mass into optical /^e light is assumed |l9i|20|| , 

Po^, = 1.4 X 10^ (^^^ Me Mpc-^ 

Here e is the overall efficiency of transforming accreted rest mass energy into 
radiation and /b is the fraction emitted in the B-band. Similarly we can 
estimate the black hole mass density which results from the emission of hard 



X-rays |21 



Px-„. = 3.8 X 10^ Mq Mpc-^ 

where we have assumed a total hard X-ray flux of 140keVs~^ cm~^ sr"'^. 
The sources producing the hard X-ray background are generally assumed 
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to be a class of AGN different to optically bright QSOs that has not yet 
been identified. We adopt an effective emission redshift Zem = f -5 for these 
unidentified sources. 

Part of the fR-background should also be produced by AGN, although it 
has been argued that their contribution should not exceed 30 percent |^ . If 
30 percent of the IR background were indeed emitted by AGN |2^, then 

ft, = 7.5 X 10^ (^^^ Mq Mpc'^ 

where we have assumed a total IR flux of I5nWs~^ cm^^ sr~^ and again 
Zem = 1.5. There has been some debate if the mass density inferred from the 
optical emission is large enough to explain the mass density in remnant black 
holes alone. Even for the low value of black hole to bulge mass ratio of 0.2 
percent adopted here, there is still a discrepancy of a factor of about five. 
This suggests (i) that there is either a significant contribution to the black 
hole mass density by accretion other than that traced by optical bright QSOs 
or (ii) that the efficiency for producing optical light during the accretion is 
lower than usually assumed or (iii) that the black hole to bulge mass ratio is 
lower than 0.2 percenct ||2j,^. The possible additional accretion may well 
explain the hard X-ray backgound and part of the infrared background. It 
may or may not trace the evolution of optical bright quasars. 

2.2 Possible assembly histories of supermassive black holes 

There is still a rather wide range of possible assembly histories for the super- 
massive black holes in nearby galaxies (see Haehnelt, Natarajan & Rees for a 
more detailed discussion |2^).In Fig. 1 we sketch three possible options out 
of this range. In the left panel most of the mass is assembled in supermassive 
black holes during the epoch of optically bright QSOs around z ~ 2.5. In 
the middle panel only 20 % of the mass is assembled in supermassive black 




Fig. 1. A sketch of possible assembly histories of supermassive black holes in nearby 
galaxies. Plotted is the overall black hole mass density relative to its present-day 
value. For furter description see text. 
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holes during the epoch of opticaUy bright QSO. The rest is accreted at low 
redshift, possibly in the form of hot gas in an advection dominated accretion 
flow. In the right panel most of the mass is assembled into small supermassive 
black holes at very high redshift. Present-day supermassive black holes form 
predominatly by merging of smaller black holes. Accretion of gas during the 
epoch of bright QSOs or at an later epoch does not change the mass density 
much. 

2.3 Clues for the formation mechanism of the typical 
supermassive black hole 

A variety of physical mechanisms for assembling mass into supermassive black 
holes have been suggested (see Rees 1984 |^ and Rees these proceedings for 
a review), 

• the dynamical evolution of a dense cluster of stellar objects, 

• the build-up of a supermassive black hole by merging of smaller black 
holes, 

• the viscous evolution and/or merger-driven collapse of a self- gravitating 
gaseous object. 

All of these processes certainly occur and can lead to the formation of su- 
permassive black holes. The last option seems, however, the most attractive 
way of explaining the observed black hole bulge mass relation. The main rea- 
son is the high formation efficiency inferred from the large black hole to bulge 
mass ratios. It is hard to see how as much as one percent of all available cold 
gas could end up in a supermassive black hole of 10^ Mq if an intermediate 
state of a dense stellar cluster is involved. Initially relaxation times in such 
a cluster would be long and a considerable fraction of stars would evaporate 
before the cluster becomes dense enough to evolve rapidly |^ . It is also prob- 
lematic to build up supermassive black holes predominantly by merging of 
smaller black holes that have formed well before the epoch of optically bright 
QSOs. This would require large black hole formation efficiencies in shallow 
potential wells with Vc < lOOkms"^. It seems more plausible that black holes 
should form less efficiently in smaller potential wells due to the feedback of 
the energy relased both by accretion onto the (forming) supermassive black 
hole and due to supernovae. It is also unclear whether supermassive black 
holes in galaxies merge efficiently or whether sling-shot ejection plays a role 
[ p8[|30t . We therefore consider the assembly history in the right panel of Fig. 
1 to be rather improbable. 

In the next section we will discuss a specific model for the evolution of 
galaxies and their supermassive black holes within a hierarchical cosmogony. 
This model assumes that the optical bright QSOs do trace the accretion 
history of supermassive black holes well (as in the left panel of Fig.l), but it 
could easily be altered to accomodate other accretion modes. 
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Fig. 2. The cosmological mass density in cold gas in galaxies as a function of 
redshift. The data is taken from Storrie-Lombardi [ pif . 

3 Modelling the assembly of supermassive black holes 
and their host galaxies 

3.1 Merging galaxies, starbursts and AGN 

In CDM-like cosmogonies, galaxies build up by hierarchical merging. The 
formation and evolution of galaxies in such cosmologies has been studied ex- 
tensively using Monte-Carlo realizations of the hierarchical build-up of galax- 
ies which include simple prescriptions to describe gas cooling, star forma- 
tion, supernova feedback and merging rates of galaxies. These models repro- 
duce many observed properties of galaxies both at low and at high redshifts 
[ pTl|l^ , [T^ , p^ . In the models, the quiescent accretion of gas from the halo re- 
sults in the formation of a disk. If two galaxies of comparable mass merge, a 
spheroid is formed and the remaining gas undergoes a starburst. We assume 
here that such major mergers are also responsible for the growth and fuelling 
of black holes in galactic nuclei. If two galaxies of comparable mass merge, 
the central black holes of the progenitors coalesce and a few percent of the 
gas in the merger remnant is accreted by the new black hole. We have made 
the following assumptions in our model: 

• The fraction of cold gas that forms stars over one dynamical time in- 
creases with decreasing redshift. 

• A fraction 0.01/[1 -f- (280/tic)^] of the cold gas in the merging galaxies 
is accreted by the black hole, where is the circular velocity of their 
combined dark matter halo. 

• The accretion timescale of the gas scales with the dynamical time, iacc — 
2.5 X 10^[0.7 + 0.3(1 + z)3]-0-5 yr. 
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(Bulge) (Bulge) 

Fig. 3. Tlie correlation between the logarithm of the mass of the central black hole 
expressed in units of W^Mq and the absolute V-band magnitude of the bulge. The 
dots are an absolute V-band magnitude limited sample of bulges in our model. The 
thick solid line is the Mv (bulge) vs. A/bh relation obtained by Magorrian et al Q 
for nearby normal galaxies. The dashed lines give an indication of the la scatter in 
the observations. The right panels is the prediction for young ellipticals. 



• A fixed fraction of the accreted rest mass energy is radiated away in the 
optical. The luminosity cannot exceed the Eddington limit. 

Our assumptions result in a strong decrease of the gas fraction in galaxies 
with redshift, from about 75 percent at z = 3 to 10 percent at z = 0. Our 
model also fits the strong decrease of the overall density of cool gas in the 
universe as inferred from the incidence rate of damped Lya absorbers (Fig 
2.). Because of this change in gas fraction with redshift, the gas fraction in 
major mergers that produce bulges is systematically higher for fainter bulges, 
which form on average at higher redshift. This might explain the systematic 
differences in the the slope of the stellar density distribution in the cores of 
high-luminosity and low- luminosity ellipticals po| . 



3.2 Remnant black holes in nearby galaxies 

Fig. 3 shows scatterplots of black hole mass versus bulge luminosity drawn 
from absolute magnitude-limited catalogues of bulges produced from our 
models. The thick solid line shows the relation derived by Magorrian et al 
and the dashed lines show the la scatter of their observational data around 
this relation. Both the slope and the scatter predicted by our models agree 
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Fig. 4. The evolution of the space density of quasars with Mb < —24. The data 
points are a compilation by Grazian et al. jssj 

reasonably well with the observed relation. The normalization is set by the 
assumed fraction of the cool gas in merger remnants which is accreted onto 
the black hole. Up to one percent of the cold gas has to find its way into 
the central supermassive black hole to reproduce a present-day black hole to 
bulge mass ratio of 0.2 percent. It is interesting to note that a central ob- 
ject with a mass fraction of about two percent can stabilize a bar instability 
in a surrounding gas disc which may be responsible for driving the gas to 
the centre pl| , ^2|j3^ . One observationally-testable prediction of our model is 
that elliptical galaxies that formed recently should harbour black holes with 
smaller masses than the spheroid population as a whole. This is illustrated 
in the second panel of Fig. 3, where we show the relation between bulge 
luminosity and black hole mass for young ellipticals. 

3.3 Evolution of optically bright QSOs 

One of the striking features of the QSO population is their rapid evolution. 
The observed rapid decline of the space density at low redshift is not trivial to 
understand |34[ . An important clue is probably the similar rapid drop of the 
overall amount of cool gas in the Universe inferred from the rate of incidence 
of damped Lya systems. The solid curve in Fig. 4 shows the model prediction 
for the evolution of the QSOs with Mb < —24 compared to observational 
data compiled by Grazian et al. |3^. The agreement is reasonably good. 
In our model the strong decrease in quasar activity at low redshift results 
from a combination of three factors i) a decrease in the merging rates of 
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Fig. 5. Host galaxy versus quasar absolute magnitudes at a series of redshifts. The 
dashed line shows the present value of L, for galaxies. The triangular box in the 
top-left panel shows the region spanned by the data set compiled by McLeod, Rieke 
& Storrie-Lombardi 

intermediate mass galaxies at late times, ii) a decrease in the gas fraction of 
galaxies with decreasing redshift iii) the assumption that black holes accrete 
gas more slowly at late times. 

3.4 QSO host galaxies 

In Fig. 5 we show scatterplots of host galaxy luminosity versus quasar lu- 
minosity at a series of different redshifts. For reference, the horizontal line 
in each plot shows present-day value of for galaxies. At low redshift. 
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quasars with magnitudes brighter than Mb = — 23 reside mostly in galaxies 
more luminous than L*. Our results at low redshift agree remarkably well 
with a compilation of ground-based and HST observations of quasar hosts 
by Mcleod, Rieke & Storrie-Lombardi The triangle in Fig. 5 marks the 
region spanned by their observational data points. At high redshifts, our 
models predict that the quasars should be found in progressively less lumi- 
nous host galaxies. This is not surprising because in hierarchical models, the 
massive galaxies that host luminous quasars at the present epoch are pre- 
dicted to have assembled recently [Q. The luminosities of quasars hosted 
by galaxies at different epochs depends, however, on the redshift scaling of 
tacc- Recently there have been a number of detection of hosts of high redshift 
QSOs 37,3^,^. Typically these seem to have ^ L, luminosity, suggesting 
that our assumed ta.cc and its scaling with redshift is indeed correct. 



4 Conclusions 

There is agreement to within a factor of a few between the black hole mass 
density infered from black holes in nearby galaxies and that inferred from 
the radiation emitted by optically bright QSOs. This is consistent with the 
possibility that optical bright quasars trace the assembly history of supermas- 
sive holes well, but significant accretion in a different accretion mode with a 
different redshift evolution is also viable. 

The large black hole to bulge mass ratio in nearby galaxies argues for 
a formation mechanism that avoids the intermediate step of a dense stellar 
cluster. A scenario in which supermassive black holes are assembled by merg- 
ers of smaller black holes which formed well before the epoch of optically 
bright QSOs would require high formation efhciencies (about 10 %) in shal- 
low potential wells. The most plausible mechanism by which the mass in a 
typical present-day supermassive black hole is assembled, is the collapse and 
accretion of cold gas plus some additional accretion of hot gas and merging 
of black holes at late times. 

It is possible to built a unified model for the evolution of galaxies their 
central black holes and AGN activity by assuming that black holes and bulges 
of galaxies form together during the frequent (major) mergers predicted by 
hierarchical cosmogonies. Such a model can reproduce the observed rapid 
evolution of the space density of bright QSO with redshift, the mass distri- 
bution of remnant black holes in nearby galaxies and the luminosity of QSO 
host galaxies. 

Interesting implications of our model are the following. The typical du- 
ration of the optically bright QSO phase should be 10'' yr. Young ellipticals 
should harbour black holes with smaller masses than the spheroid popula- 
tion as a whole. QSO hosts are typically brighter than at low redshift and 
should become fainter with increasing redshift. Important for the rapid de- 
cline of the space density of bright QSOs and the cosmological density of cold 
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gas, is that the gas fraction in galaxies decreases with decreasing redshift. As 
a consequence, fainter ehipticals have formed in more gas-rich mergers than 
bright eUipticals. Supermassive binaries and merging of supermassive bina- 
ries should occur frequently in hierachical cosmogonies. The latter is good 
news for space-borne gravitational wave experiments like LISA |^,^ . 

We acknowledge helpful discussions with Andrea Cattaneo, Stefano Cris- 
tiani, David Merrit, Prija Natarajan, Joel Primack, Martin Rees, Hans- Walter 
Rix and Simon White. 
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